function y = xyvalcs(p1,p3,psis,NN,NLS,wfr,wfi,wfr1,time1,time2)at=p1(1:NLS);bt=(wfr*real(p1(NLS+1:2*NLS))+1i*wfi*imag(p1(NLS+1:2*NLS)));ct=wfr1*p3(1:NLS);bo=at;y=zeros(1,NN);y(1)=at*psis(:,1);for nr = 2:NN    gc = exp((real(bt).*abs(time2(nr)))+(1i*imag(bt).*time2(nr)));    %  gc = exp(bt.*abs(time2(nr)));    co = exp(ct.*time1(nr));    y(nr) = ((bo.*gc).*co)*psis(:,nr);end